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Abstract 

The velocity potential in the Kelvin ship-wave source can be partly expressed in 
terms of space derivatives of the single integral 

/ OO 

exp [—ipcosh(2u — ia)] cos(a;coshu) du, 

-OO 

where {x, p, a) are cylindrical polar coordinates with origin based at the source and 
— ^Tr<a< ^TT. An asymptotic expansion of F{x, p, a) when x and p are small, but 
such that M = x^ jiAp') is large, was given using a non-rigorous approach by Bessho 
in 1964 as a sum involving products of Bessel functions. This expansion, together 
with an additional integral term, was subsequently proved by Ursell in 1988. 

Our aim here is to present an alternative asymptotic procedure for the case of 
large M. The resulting expansion consists of three distinct parts: a convergent 
sum involving the Struve functions, an asymptotic series and an exponentially small 
saddle-point contribution. Numerical computations are carried out to verify the 
accuracy of our expansion. 
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1. Introduction 


The velocity potential in the Kelvin ship-wave source can be partly expressed in terms 
of space derivatives of the single integral 

/ OO 

exp [—\z cosh 2u + ^iy sinh 2u] cos(x cosh u) du 

-OO 


exp 


[—i/3cosh(2u — ia)] cos(x coshu) du = F{x, p, a) 


( 1 . 1 ) 

where the origin of cylindrical polar {p, a, x) coordinates is at the source, x > 0 is the non- 
dimensional horizontal coordinate along the track, y = p sin a is the horizontal coordinate 
perpendicular to the track and z = p cos a is the vertical coordinate (increasing with 
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depth). The polar angle a is zero directly below the origin and satisfies — ^vr < a < 
although it is sufficient to restrict attention to 0 < a < since F{x,p,a) is an even 
function of a. 

The integral F{x, p, a) is difficult to evaluate when x and p are small. In 1964, Bessho 
[1] gave the expansion 

OO 

F{x,p,a) = Ko{^p)Jo{x) + 2 ^ cosma Km{lp)J 2 mix), 

m=l 

which he obtained by expanding the cosine factor as a series of Bessel functions. Here J 
and K denote the usual Bessel functions, and from the well-known large-order behaviour 


J2m{x) ~ 




ep \ 
4m / 


(m ^ -|-oo) 


it follows that the late terms in the above sum are controlled by {ex"^/Ap)'^/. The 
expansion is therefore absolutely convergent, but in the case where x and p are small 
with x‘^/{Ap) assumed to be large, the convergence is slow and the accuracy is reduced 
by the accumulation of cancellation errors. 

An expansion suitable for M = x^/{Ap) 1 was established rigorously by Ursell [7] 

in the form 

F{x,p,a) = -7r(lo{lp)Yo{x) + 2 ^ cosma/m(^/o)> 2 m(a:) 

^ m<M 


/ —oo-\-^i\a.\—7vi 

exp [—i/9cosh(2u — i\a\) + ixcoshu] du + 0{e~^), (1.2) 

-oo+ii|o| 

where I and Y are the usual Bessel functions. The series extended over infinite values of 
m (without the integral and order term) had been obtained in a non-rigorous fashion in 
[1]. From the large-m behaviour 


Im ioP) 


y/2TTm \AmJ 


Y2m{x) ~ 




(m ^ +oo), 


it is easily seen that the terms in the series in (1.2) are comparable to an expansion in 
inverse powers of M. The additional term represented by the integral in (1.2) can be 
estimated by the method of steepest descents and is given by [7, Eq. (2.26)] 

(7r/M)5e-(^-5^)“"“ sin[(M + ip)sina+ia], (1.3) 

which is exponentially small except when a is near ^tt. An algorithm for a more precise 
estimation of this term is given in [4]. 

Our aim in this note is to present a different asymptotic procedure for F{x.p, a) when 
X and p small, but with the parameter M 1. We deal first with the expansion in the 
midplane (that is when a = 0) and then extend the procedure for 0 < a < ^vr. The 
expansion we obtain consists of three parts; the first is a convergent series involving 
Struve functions that is rapidly convergent for small values of x and p; the second is 
an asymptotic series and the third is an exponentially small saddle-point contribution. 
This latter contribution agrees with (1.3) when a is bounded away from zero. Numerical 
calculations reveal that our expansion and that in (1.2) produce comparable accuracy. 
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2. Asymptotic evaluation when a = 0 


We consider the case where x and p are small but x‘^l{4:p) is assumed to be large. When 
a = 0, we have from (1.1) 




/ OO , , fOO p — lVlT 

cos(xcoshtt) du = 2e2^ / - j 

-oo Jp Y ' 

where we have set coshtt = xTf{2p) and defined 

2p 


cos (2Mr) 


t 2 _ p2 


dr, (2.1) 


M = 


4p 


p = — = 


x 

i/if' 


The integrals 


foo M(r^^2jr) 


dr 


^ -j ... (2.2) 

Jp V T - P 

have saddle points at r = ±i, with the steepest descent paths through these points being 
the doubly infinite lines parallel to the real r-axis. The parts of these paths situated in 
3?(r) > 0 are labelled Li and L 2 in Fig. 1. The square root in the denominator makes it 
necessary to introduce a branch cut the r-plane along [—p,p]. 

Considering the integral with the upper sign in (2.2), we deform the integration path 
Ip, 00 ) into a semi-circular path round the branch point t = p, the upper side of the 
branch cut to the origin followed by the positive imaginary axis 0 < 9 (r) < 1 and then 
out to infinity along the steepest descent path (which we label Li) through t = i. This 
corresponds to the paths shown in the upper half of Fig. 1. If the phase of (r^ — is 

zero on the positive axis on the right of r = p, then the phase is on the upper and 
lower sides of the branch cut, and on the positive and negative halves of the imaginary 
axis, respectively. Then we find, on letting the radius of the indentation around t = p 
shrink to zero, 


^00 „—M(T^— 2 iT) 


_ p2 


■ dr = i 


fP ^-M{T^-2iT) 

0 \/p2 _ -^-2 


dr + 


fl M(u^—2u) 


'0 \Jv?‘ + p'^ 


du — i 


AL(t^—22t) 


Ll \/p2 _ t2 


■ dr 


Similar considerations applied to the integral with the lower sign leads to the above 
relation with i replaced by —i and the path Li by L 2 . Then we find 


■2^F(x,p,0) = -2 / 
Jo 


P e sin (2Mr) 

\/p2 _ j2 


dr + 2 


fl ^M{v?—2u) 


0 +p^ 


du + Is 


where Ig denotes the saddle-point contribution 

^-M{T^+2iT) 


Is = i 


IL 2 


- 


dr — i 


g—AL (t^— 2iT) 


/Li 


Vp"^ - 


dr. 


(2.3) 


(2.4) 


We now consider in turn each of the integral contributions in (2.3). For this we shall 
require the evaluation^ 

rg-MF^2fc+i(^2 _ + i;k + p + l-, -Mp^) (2.5) 

Jo 2T{k + pL + l) 

^This result can be obtained by using the series expansion of the exponential and evaluating the 
resulting integrals in terms of the Beta function. The result given in [3, 3.478(3)] incorrectly excludes 
this case. 
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Figure 1: The deformation of the integration path [p, oo) in the r-plane through the saddles at 
±i. There is a branch cut along [—p,p]. 


for /i > 0 and k = 0,1,2 ,where iFi denotes the confluent hypergeometric function. 


2.1 The branch-cut contribution 

We first consider the branch-cut contribution 



sin (2Mt) 

^p2 _ ^2 


dr. 


Expanding the sine function followed by reversal of the order of summation and integra¬ 
tion, we find 


h 


E 

k=0 


— fP .j-2k+l 


(2A:-|-1)! Jo 


dr 





(-)^(x/2)2fc+i 

r2(A: + |) 


■j- iFi{k + 1; A; + |; —p) 
2 ) 

Fi(^;A: + |;p) 


by (2.5) with /r = i, Mp^ = p and 2Mp = x. Here we have employed Hummer’s 
transformation for the confluent hyper geometric function and also the duplication formula 
for the gamma function. 

If the iFi function is expressed in its series form, we then obtain 


h 


TTC 


-n 00 (1 




(-)‘( U .)“+1 


r=0 


r(/^ +1 )r(fe + r + f) 


Tre 


-o 00 (1 




r=0 


rp I 


( 2 . 6 ) 


where {a)r = r(a -|- r)/r(a) is Pochhammer’s symbol and TLu{x) is the ‘scaled’ Struve 
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function defined by [5, p. 288] 


'Hu{x) 


OO 


(ix) 

k=0 


r(fe + |)r(fe + p + |)’ 


(2.7) 


The series on the right-hand side of (2.6) is absolutely convergent, since the behaviour 
of the scaled Struve function as r —)■ oo is given by Tir (x) ~ xe'’/{^^^/2r''^^); see [5, 
Eq. (11.6.5)]. 


2.2 The imaginary axis contribution 

We now consider the second integral in (2.3), namely 


h 


nl pM{v?—2u) 


/o \Ju^ + 


du. 


This represents the contribution from the origin down the steepest descent path (relative 
to the origin) to the saddles at r = and is a Laplace integral with a linear endpoint at 
u = 0. However, it is not possible to employ the standard asymptotic estimate for such 
an integral as described in, for example, [5, 2.4(iii)], [6, p. 14], since we have M —)• oo 
and p = x/{2M) —)■ 0. With the change of variable ^ = u/p, this becomes 


h 




/o yr+F 


dt 


1 



2M 

X 


We now write as a finite series of n (> 1) terms together with a remainder 
expressed in Lagrange’s form [2, p. 328], [9, p. 96]: 




2\k 


= E 

k=0 


k\ 


+ 




2\n 


jpe 


n\ 


(0 < 0 < 1 ). 


This then enables us to write 



^2kp-xi 

7^ 


d^ + Rn, 


( 2 . 8 ) 


(2.9) 


where 


Rn = Y , 

nl Jo 


pTi .^0 ^‘2.np-xi+epi'^ 




di. 


The main contribution to the finite sum in (2.9) may be written in the form 


( 2 . 10 ) 


Q - V 
- 9 Z- 22fcA:! 


7C'fc(x), Cfc(x) := 

' tt Jo 


CO 


di. 


To evaluate the integrals C'fc(x), we observe that for integer m > 0 [5, Eq. (11.5.2)] 


( 2 . 11 ) 


E 

r=0 


r tS? K„(x), 


Kpiz) = Hp{z) - Yp{z), 
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where K.^{x) and ^y{x) are Struve functions and Yy{z) the usual Bessel function; see [5, 
Eq. (11.2.5)]. Then, with /Cm(x) = we obtain the recurrence relation 

C^{x) = 2™(i)^/C™(x) - E i ^)a(x) (m > 1), (2.12) 

r=0 ^ ^ ^ 


from which we find the first few C'fc(x) given by 

C'o(x) = /Co(x), C'i(x) = /Ci(x) - x^/Co(x), 

(72 (x) = 3/C2 (x) — 2x^/Ci(x) + x^/Co(x), 

(73(x) = 15/C3(x) — 9x^/C2(x) + 3x^/Ci(x) — x®/Co(x), 

(74 (x) = 105/C4 (x) — 60x^/C3(x) + 18x^/C2(x) — 4x®/Ci(x) + x®/Co(x) . 

In Appendix A it is established that the contribution resulting from the integral over 
[.^ 0 )Co) in (2.9) is 0{e~^), when^ n < M, and that Rn = Then, from (2.11) 

we have the result 


h = l y. 


n—l 


M-^ 


CJx) + 0(M-^] 


(2.13) 


as M —)■ oo. 


2.3 The saddle-point contribution 

The final term to consider in (2.3) is the saddle-point contribution R defined in (2.4). 
Recalling that the integration paths Li and L 2 are r = tt±i, 0<u< 00 , we have 


/o = ie 


-M 


e-^^^f{u)du, 


fiu) = 


— {u — iY a/p^ — {u + iY 


Then, since f{u) = 2m/(l -|-p^)^/^ -|- O(u^) as u —>■ 0, we obtain 


2e~^ T^ 
= (l+p2)3/2 I 


^-Mu^ 


{u -\- (7(tt^)} du = 


,-M 


M(1 -|-p2)3/2 


{1 + 0{M-^)} (2.14) 


as M —)■ 00 . 


From (2.3), (2.6), (2.13) and (2.14), we obtain our desired result when a = 0 given 
in the following theorem. 

Theorem 1. Let x and p he positive and bounded (to fix ideas suppose 0 < x < 1 and 
0<p< 1) and let a = 0. Then, when M = x^/(4p) is large we have 


00 /1' 


F(x,p,0) =-vre z/’^ /R^(x)-h 7re2/’| ^ 


r=0 


n-1 ^-k 


22kk\ 


Ckix) -\- 0{M 


'k=0 

-—M\ I n/r—l „—M\ 


+0(e"^'") + OiM-^e-^^^) (2.15) 

where n < M is a fixed positive integer, Tir{x) is a scaled Struve function and the 
coefficients Ck(x) are defined by the recurrence in (2.11). 


^If n is finite, this estimate can be replaced by 0{M"e see (A.4). 
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3. Asymptotic evaluation when 0 < a < ^vr 

When 0 < a < ^vr, we can make the substitution u ^ u + ^ia in (1.1), appealing to 
Cauchy’s theorem to place the integration path back onto (— 00 , 00 ), to find 

/ OO ^ 

g- 2 Pcosh 2 « cos(xccoshtt -|- isxsinhu) du, 

-OO 

where, for brevity in what follows, we have set 

c = cos ^a, s = sin ^a. 

With coshtt = xt/{2p), some routine algebra then shows that 

e~2^F{x, p, a) = J^{M, a) + J-{M, a). 


where 


J±{M, a) := 


,M( 


-F±2icT) cosh(2MsVr^ -p2) 


- 


dr. 




(3.1) 


The integral J+(M, a) contains two components associated with the exponential fac¬ 
tors exp [Mi/) (r)], where 


i/)(r) = —-|- 2icT ± 2s\J 

we shall refer to these as the ‘positive’ and ‘negative’ components of J+{M,a), respec¬ 
tively. The components are associated with saddle points given by 


r, = ic±s^-^^^. (3.2) 

Remembering that p = x/{2M), we see that for large M each component of J^{M,a) 
has a saddle situated at 


rt (p ^ 0). (3.3) 

Some routine algebra, using (3.2) and (3.3), shows that 

^ 2icTt + ~ _ t -±2 ^ ^ ^p2^^±ia/2 

Ts - ic 

~ ± (p^o). (3.4) 

Typical paths of steepest descent through these two saddles are shown in Fig. 2(a), 
where that through the saddle t~ is shown dashed; only the saddle will be of relevance 
in our considerations. Both these paths, which we label and respectively, intersect 
the imaginary r-axis at the common point r = mo, where 

’^{'4){Tf)) = ±.2s^Uq -|-p2 ~ ±(sina-|- ^p^sc); 

that is 

no — c(l — tan^ ^a) (0 < a < ^tt, p ^ 0). (3.5) 
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Figure 2; (a) The steepest descent and ascent paths through the saddles for the two 
components of J+{M,a) and the point iug when P = \ and a = 27r/5. The steepest paths 
through r+ in 5R(r) > 0 when (b) a = 7r/4 and (c) a = 7r/2. There is a branch cut along [—p,p]. 


The integral J-{M, a) is associated with a conjugate set of saddles in the lower half-plane 
and conjugate steepest descent paths and When a = 0 the steepest descent paths 
become the horizontal lines Li and L 2 emanating from the points r = shown in Fig. 1. 

We now deform the integration path as described in Section 2, recalling that the 
phase of (r^ — is on the upper and lower sides of the branch cut [—p,p]- 

The integral J+(M, a) is taken along the upper side of the cut to the origin, followed 
by part of the positive imaginary axis [0, iuo] and then each component is taken out to 
infinity along its associated steepest descent path Cf and Cj". The positive component 
of J+(M, a) passes over the saddle at r+, whereas the negative component passes over 
part of the steepest descent path situated in 5R(r) > 0 that emanates from the saddle t~ 
in 5R(t) < 0. The integral J-{M, a) is taken along a similar path in the lower half-plane 
passing out to infinity along the paths and C 2 ■ Then we find 


■2^F(x,/9,a) = -2 / 

Jo 


' -Mt'^ ■ \ COs(2Ms\/p2 _ t-2^ 

e sm(2Mcr)- J ^ - - dr 


Vp'^ - 




+2 


^M{u^-2cu) COs{2Ms^Jv? +p‘^) 
0 yjv? +P^ 


du + Is 


(3.6) 


where 

j oc- /■ ^-M{T^-2icT) cosh(2Msv^^^ 

J C \/p^ — 

and by the path C it is understood that this refers to the steepest descent paths Cf^ 
situated in 3?(r) > 0 for the positive and negative components of the integral, respectively. 
When a = 0, the expression on the right-hand side of (3.6) agrees with (2.3) and (2.4). 
In the hrst integral in (3.6) we expand the trigonometric functions to obtain 


h = 


V V { 2 Msf^ 


) m=0 
— o 00 

7re ^ ^ 


(2m)! 


^\2m np 

^)! Jo 


00 00 / 






upon use of (2.5) with p = m+^ and application of the duplication formula for the gamma 
function and Rummer’s transformation. Then substitution of the series expansion of the 
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hypergeometric function yields the convergent expansion 


h 


(3.7) 


where ^^{x) is the scaled Struve function defined in (2.7). When a = 0, only the term 
corresponding to m = 0 makes a contribution and we recover the evaluation given in 
( 2 . 6 ). 

With the change of variable ^ = u/p, the second integral in (3.6) can be written as 


h 

where the value of mq as p —)■ 0 is given in (3.5). We follow the same procedure employed 
in Section 2.2 to find (see Appendix A for details) 


pe-xci cos(sxyi + g^) 

Jo vTTF 


Co = 


2Muq 


h = £ i r ^ ^ o(M-) 

fc=o -^0 VI + 

^ 71-—1 l\/f _ k 


when n < McP, where 


^ / \ 2 2fe /■°° ^2k -xcf cos(sxa/ 1 -k C^) 

Cfc(x,a) = -x^M - /rr^ 

71 Jo VI +C 


di 


2 .oj, cos(sVx2 + V) 


t^^e 


vr Jo 


Va7^ + V 


dt; 


(3.8) 


(3.9) 

(3.10) 


compare (2.11) in the case a = 0. 

Finally, we determine the contribution Ig- The positive component of the integral Ig 
passes along the steepest descent path through the saddle r^, where for a bounded away 
from zero we obtain by application of the saddle-point method the estimate 

( 27r y ^Mi;(T+)+hia 

V-MV"(rs+); 

where ~ —2. From (3.4), we have 

V(t^) — —-|- (p —)■ 0). 

The contribution from the negative component is more recessive (being of 0(M“^e“^‘^^)), 
since it involves a portion of the steepest descent path not including the saddle point 
T~. Hence we find the estimate for large M 

e^PIg ~ (7r/M)ie-("-^^)“"“ sin [(M + ip) sina + ia] (3.11) 

when 0 < a < ^tt, which agrees with that obtained by Ursell stated in (1.3). This 
estimate is exponentially small except when a is near ivr. It is also important to stress 
that the above estimate is not valid as a —)■ 0; see (2.14) when a = 0. 

From (3.6), (3.7), (3.8) and (3.11) we obtain our desired result when 0 < a < ^vr 
given in the following theorem. 
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0.0 0.5 1.0 1.5 2.0 


Figure 3: Behaviour of the first few coefficients Ck{x,a) as a function of x when a = tt/G. 


Theorem 2. Let x and p be positive and hounded and let 0 < a < ^vr, with s = sin ^a, 
c = cos ^a. Then, when M = l(A,p) is large, we have 


F{x,p,a) 


—Tre 


oo r 

-hpy^ 

r=0 ’ m=0 


^ (-r(m + i) 


ml 


{ysyLLm+rixc) 


1 ^ 


(3.12) 


where n (< Mc^) is a fixed positive integer, ^^{x) is the scaled Struve function defined 
in (2.7) and the coefficients Ck{x,a) are defined in (3.10). 


Some examples that demonstrate the precision of the asymptotic formulas in Theo¬ 
rems 1 and 2 will now be given. In Fig. 3 we show the behaviour of the first few coefficients 
Ck{x,a) as a function of x for a given a. Denoting the convergent sum involving Struve 
functions in (3.12) by S'! for brevity, then we define 


F := F{x, p, a) + vre '^^Si — 


n—1 


M-^ 


Table 1 shows the absolute error in the computation of F from the expansions in (2.15) 
and (3.12) for different values of a together with the truncation index n employed (satisfy¬ 
ing n < Mc^). The convergent sum Si was computed with sufficient terms commensurate 
with the overall accuracy of the computation. It is found that these errors are compara¬ 
ble to those resulting from optimal truncation of the Bessho expansion in (1.2) (without 
the integral term). 


Appendix A: Derivation of bounds on the contribution from the imaginary 

axis for 0 < a < y 


The contribution to F(x, p, a) arising from the imaginary r-axis over the interval [0, mo] 
is expressed in the form 



cos(sx\/r+^) 


df, + Rni 


(A.l) 


vTTF 
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Table 1: Values of the absolute error in the computation of T and the truncation index n. 



X = 0.40, p = 

0.005 

X = 1.00, p = 0.020 


(M = 8) 


(M = 12.5) 


a/i: 

|R| 

n 

\F\ 

n 

0 

6.368 X 10"® 

8 

2.613 X 10-'^ 

12 

0.10 

3.146 X 10"® 

5 

1.998 X 10-® 

12 

0.20 

3.146 X 10"® 

6 

1.899 X 10-^ 

11 

0.25 

4.687 X 10-3 

5 

1.428 X 10-3 

9 

0.30 

2.976 X 10-3 

3 

2.890 X 10-^ 

9 

0.40 

4.326 X 10-2 

1 

7.928 X 10-^ 

8 


where 


-Ru 


^ -xci+eoe COs(sXyi + 

n\ Jo ^ vT+F 


0 < 0 < 1 


and ^0 = 2Muojx. This results from expansion of as a finite series of n terms 
together with a remainder in Lagrange’s form; see (2.8). The quantity uq is the intercept 
on the imaginary axis of the steepest descent path through the saddle labelled From 
(3.5) we have uq = c(l — ip^tan^ ^a) as p —)• 0, so that uq < c. When a = 0, we have 
uq = 1, ^0 = ‘2'Mc/x and the above expressions reduce to those given in (2.9) and (2.10). 

The remainder term R„ satisfies 


\Rn\ < ^ / = ^ ^2n-l^-xci/2 . 

n\ Jo n\ Jo 

where g{^) = exp [—+ 6p^‘^]. It is easily shown that g{^) is monotonically decreasing 
on [0,^o] when 0 < 0 < 0O) = c/{2uo), and that there is a single minimum at ^ = 

Me/{Ox) when 6q < 9 < 1. Hence it follows that g{^) < 1 on the interval [0,^o] since, 
with uo < c, 

^(^o) = < 1. 

Then we have the bound 




n! Ji 


2n-l^-xcil2^^ 


< 


M-^T{2n) 


n\ 


(A.2) 


upon replacing the upper limit of integration by oo. 

The tail T of the finite sum in (A.l) involving the integral over [^QjOo) can be esti¬ 
mated as follows: 


f nk roo p2kp—xd^ u 1 h 




fc =0 
n—1 


= E 


{Me 


.2\-k 


k=0 


22kk\ 




n—1 


f 2 Muqc 


dw=J2 

k=0 


{Me 


2\-k 


22kk\ 


T{2k, 2Muoc), 


where r(a, y) denotes the upper incomplete gamma function. We now suppose that 
n < Muoc < Mc^ and employ the simple bound (see Appendix B) for a > 0, y > 0 given 
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by r(a, x) < 2x“e ^ when 0 < a < x- This yields the bound, for n < Muqc < Mc^, 




^ 2g“2Mtioc (-MUq) ^ 2g—2Muoc+Muf^ 


k=0 

< {n<Mc^). 

If n is finite, the above bound can be replaced by 

n—1 

|r| < 2e"2^“°‘^^(Mug)'= = 

A:=0 


(A.3) 


(A.4) 


Appendix B: Derivation of a bound on r(a, x) 

The derivation of this bound follows closely that given in [8]; see also [6, p. 76]. Suppose 
that a > 0 and x > 0- Then [5, p. 174] 

roo roo 

r(a + 1, x) = / e-H'^dt = / e-“^(l + ufdu. 

Jx Jo 

When — 1 < a < 0, we therefore have 

roo 

r(a + 1, x) < / e~^^du < ix > !)• 

Jo 

When 0 < a < 1, we have (1 + u)“ < 1 + tt for u > 0, so that 

r(a + l,x) < e“^x“^^ / e““^(l + u) du < e^x“^^ f—+ 

< 2e->^x“+^ (X > 1). 

When a > 1, we have 

roo 

r(a + l,x) < e-^x“+W {e-“(l + u)}“du 

^0 

roo 

< + u)du = } 

Jo 

the first inequality holding when a < x and the second inequality when a > 1, since the 
integrand e““(l + u) < 1 on [0,oo). The resulting bound therefore holds for 1 < a < x- 
Collecting together these results and replacing a by a — 1, we therefore obtain the 
upper bound 

r(a, x) < 2e"^x“ (0 < a < x) (B.l) 


provided x > 1- 
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